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The time evolution of the fluorescence intensity emitted by well-defined ensembles of Green Flu- 
orescent Proteins has been studied by using a standard confocal microscope. In contrast with 
previous results obtained in single molecule experiments, the photo-bleaching of the ensemble is 
well described by a model based on Levy statistics. Moreover, this simple theoretical model allows 
us to obtain information about the energy-scales involved in the aging process. 
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Since the beginning of the 90's, single molecule exper- 
iments have opened the way to the study of the optical 
properties of a single emitter 0. In particular, these 
experiments have evidenced that an aging mechanism af- 
fects fluorescent objects: each marker performs only a 
finite number of absorption/spontaneous emission cycles 
before definitely reaching a permanent photo-bleached 
dark state Q. Before being irreversibly photo-bleached, 
single emitters often display a "blinking" behavior: as 
the time goes on the system undergoes transitions be- 
tween a bright state (ON state in which it stays a time 
t n) and a reversible dark state (OFF state in which it 
stays a time tqff)^- hi a single molecule experiment 
one can measure the light emitted by a single marker as 
a function of time and then characterize the blinking and 
photo-bleaching behaviors. The study of a multitude of 
single systems gives access to the probability distribution 
associated to their emission properties, for instance the 
statistical distribution of emission spectra 0, |f| or the 
probability distribution P(t on /off) associated to ton 
and toff 0- 0- I n this way it is possible to obtain 
more complete information with respect to an ensemble 
measurement. Such a strategy works well in the case of 
inorganic emitters (i.e. semiconductor quantum dots). 
In this case the observation time of a single object (lim- 
ited by photo-bleaching) is long enough to obtain a large 
number of measurements and therefore meaningful his- 
tograms. Moreover, such systems allow the implementa- 
tion of parallel measurements on some hundreds of sin- 
gle emitters (wide field microscopy) opening the way to 
study the ergodicity of such systems 0, Q ■ Several ex- 
periments show that the ton and toff are distributed 
according to broad laws lia] • Let us stress that a long 
observation time is crucial for a complete characteriza- 
tion of such broad distributions, due to the intrinsic large 
dispersion of the measured data. In the case of single or- 
ganic molecules, the observation time is strongly reduced 
by quenching phenomena. This fast photo-bleaching is 
a limiting factor for the obtention of reliable statistical 
distributions except in special cases of well controlled 
environment @. A particularly interesting case of or- 



ganic markers is that of fluorescent proteins such as the 
Green Fluorescent Protein (GFP) which is one of the 
most popular marker for in vivo imaging in cell biology 
[Io|. Single molecule room temperature studies on GFP 
show results that are compatible with an aging mecha- 
nism governed by narrow statistical laws (exponential or 
bi-exponential). With this assumption it is possible to 
determine average values such as (ton), (toff) and the 
characteristic time associated with the irreversible tran- 
sition between the fluorescent and the dark state (the 
average photo-bleaching time) Reported average 

photo-bleaching times for fluorescent proteins are of the 
order of some hundreds of millisecond, depending on ex- 
citation intensity Therefore single molecule experi- 
ments seem to evidence a fundamental difference between 
inorganic and organic emitters: the times associated to 
aging are distributed according to broad versus narrow 
laws. However, the presence of Levy statistics cannot be 
excluded if we consider the poor quality of the histograms 
obtained in the case of proteins. 

The question that we address in this letter is : Can one 
find signatures of Levy statistics in ensemble measure- 
ments in the case of fluorescent proteins? The analysis 
of the time evolution of the fluorescence intensity emitted 
by a well defined ensemble of GFP under c.-w. excita- 
tion shows that this is the case. Starting from ensemble 
measurements, we evidence that the GFP aging process 
is governed by Levy statistics. Moreover, information 
about the potential barriers involved in the aging phe- 
nomenon is extracted from the experiment, thanks to a 
theoretical model that takes into account an irreversible 
photo-conversion between a fluorescent and a permanent 
dark state. One of the main results of this study is that 
the unstable bright state has an infinite mean lifetime. 

In the experiment, conceptually very simple, we have 
recorded the time evolution of the fluorescence intensity 
F emitted by a well-defined ensemble of a GFP mutant 
(uv mutant : GFPuv |l3jl under continuous excitation. 
The sample is obtained by dispersing a recombinant GF- 
Puv solution (125 fjM in a pH=8 buffer solution) in a 
polyacrylamid matrix To avoid denaturation, the 
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protein solution is added after polymerization. The fi- 
nal concentration in the anhydrous samples is around 
4 x 10 5 molecules per /im 3 . The sample is placed on 
a slide-holder of a home-built confocal microscope illu- 
minated by an Ar-ion laser (A=457 nm, 0.01<intensity 
I <2 MW/cm 2 ). A microscope objective of numerical 
aperture 0.65 (x40) leads to a theoretical lateral resolu- 
tion of approximately 300 nm. The detection pin-hole 
has a radius r p i n of fO [im. Our ensemble is then de- 
fined by the ~200000 immobile molecules present in a 
detection volume of 0.5 [irn? . The emitted fluorescence 
F is detected by a photo-multiplier in current mode cou- 
pled to a lock-in amplifier. Let us stress that, in order to 
have reliable results, the ensemble of excited molecules 
must not be renewed during the whole observation time. 
We have tested the long-term mechanical stability of our 
setup and verified that it is of the order of the lateral 
resolution of the microscope. 
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FIG. 1: Normalized fluorescence intensity F(t)/F(t — * 0) 
(solid line) emitted by a well defined ensemble of GFPuv 
(~ 200, 000 proteins) as a function of the total illumina- 
tion time t. Experimental conditions: J=0.96 MW/cm 2 and 
T=290 K. In dotted line the power law used to fit the long 
time tail: (-^-) a with ta = 6.8 s and a = 0.73. For com- 
parison we report also a thermal aging curve obtained at 
T = 353 K (circle) in a fluorimeter together with a dotted-line 
curve which represents its best exponential fit. 

In Fig. 1 we present the observed time evolution of F 
as a function of the total illumination time t. All the en- 
sembles that we have studied under different illumination 
conditions display a time evolution of F(t) with a long- 
time tail well described by a power law (-y-) Q character- 
ized by an exponent a < 1. Let us stress that in order to 
unambiguously identify the asymptotical power law be- 
havior of F(t) it is compulsory to scan several decades in 
time with a good detector dynamics. The aging described 
by the curve F(t) appears to be an irreversible phe- 
nomenon that depends on the illumination history of the 
ensemble. To investigate this phenomenon, we switchcd- 
off the illumination at time to and kept the sample in 
the dark during a time Tdark- By switching-on again 
the illumination, F recovered its value at F(to) for Tdark 
up to 24 hours. We conclude that the unique relevant 
parameter for the description of the aging curves is the 



total illumination time t. These observations show that 
the transition between the fluorescent and the permanent 
dark state takes place in the electronic excited state of the 
protein. One may argue that the observed aging curve is 
the consequence of the thermal denaturation of the pro- 
tein induced by laser-heating ^(J . We have excluded this 
hypothesis by performing several measurements of ther- 
mal denaturation in a temperature-controlled cuvette un- 
der low intensity illumination (uv lamp in a fluorimeter) . 
A significant thermal denaturation of GFPuv is clearly 
observed for temperatures higher than 340 K, but this 
phenomenon induces a purely exponential loss of the flu- 
orescence as a function of time (see Fig. 1). 

In order to describe the experimental results, we pro- 
pose a model that postulates the presence of Levy statis- 
tics in the irreversible transition between the fluores- 
cent and the dark state of the protein. Within this 
model, photobleaching is the consequence of the ther- 
mally activated crossing of a barrier in the electronic ex- 
cited state of the protein (see scheme in Fig. 2). We 
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FIG. 2: Scheme of the relevant energy-levels and coupling 
constants used to model the aging mechanism of the GFPuv. 
\g) and |e) are the electronic ground and excited states of the 
protein, while \D) represents a permanent dark state. \e) and 
\D) are separated by a thermally-activated barrier of height 
E. r, r' and Pi = - are, respectively, the spontaneous emis- 
sion, the absorption (and stimulated emission), and the con- 
version towards the dark state rates. 

make the hypothesis that the heights of the barriers in a 

sample are distributed according to an exponential law 
i — — 

P(E) = -g-e B ° , where Eq is the characteristic height 
of the barriers. The time Tj associated to the crossing of 
the i th barrier of height Ei is then given by the Arrhe- 

nius law T\ = Toe k B T , where r is a characteristic time, 
kg the Boltzmann constant and T the temperature. By 
considering that P{E)dE = P(r)dr, it is easy to show 
that, for t » r , 

P{r)=^, (I) 

with a = 15]. With these assumptions, the asymp- 
totic fluorescence at time t is simply proportional to the 
number of proteins that have a bright state lifetime (time 
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spent by the protein performing absorption/spontaneous 
emission cycles) greater than the time t. The fluorescence 
F is thus proportional to the survival probability of the 
molecules in the fluorescent state: F(t) oc P(r)dr. 
At long times we expect that F(t) oc (tP-)™, which exactly 
corresponds to the experimental observations. Within 
this model, the fact that the observed a coefficient is 
lower than 1 has a fundamental consequence: the exis- 
tence of an unstable bright state with infinite mean life- 
time. The bright state is unstable in that, at long times, 
each molecule ends up in a permanent dark state with a 
probability tending to one. The mean lifetime is infinite 
in that the ensemble average of the time r (time of jump 
to the dark state or photobleaching time) is infinite. 

Let us develop now the model more in detail. The three 
relevant states (scheme in Fig. 2) are coupled by the ab- 
sorption and the stimulated emission (rate V oc I), the 
spontaneous emission (rate T), and the bleaching process 
which ends in the dark state \D) (rate T\ = i). By fix- 
ing the excitation intensity / ( V) and the rate Ti, we 
can solve a system of coupled differential rate equations 
that gives us the time-dependent population ir e (Ti, I,t) 
in the |e) state. As long as we are not interested in the 
short-time behavior (fast equilibration of electronic pop- 
ulations) we obtain, for T\ << (r + 2r'), the simplified 
expression: 



7r e (r,t) 



r' 



r + 2r 



(2) 



In order to compare the model with the experiments, 
two fundamental ingredients have to be added to it: the 
spatial distribution of the excitation intensity (typically 
a gaussian laser spot with characteristic dimension ro 
/ = I{r) = 2o e ) an d the stochasticity associated to 
the rate T\. 

The first aspect has been discussed by A. J. Berglund 
in a theoretical paper 0] that addresses the case of the 
samples illuminated by intensity profiles which span an 
infinite object plane. In this case, a gaussian laser spot 
and a simple model for the photobleaching (e.g. Ti rate 
constant) lead to a time evolution of F with a long time 
tail that decays according to the power law F(t) oc | [l7| . 
At first sight, this result seems to apply to our experimen- 
tal observations. However in our experiment (and, more 
generally, in the experiments performed in a confocal ge- 
ometry) the detection volume is limited by the pin-hole 
aperture: the pertinent intensity distribution is therefore 
a truncated gaussian. The consequence of this trunca- 
tion is that, for T\ distributed according to a narrow dis- 
tribution law, F(t) displays an asymptotic exponential 
behavior. We have checked by numerical integration of 
equation © that, for realistic sizes of detection pin-holes 
(i.e. r /2 < r pin < 2r ), the spatial distribution of the 
excitation intensity can only influence the details of the 
curve F(t) at intermediate times but not its asymptotic 
behavior. 



The properties of the statistical distribution of Ti seem 
therefore to be the key element to understand the ob- 
served behavior. In order to introduce this stochastic- 
ity, we have to consider the ensemble average of ©: 
(ir e (t)) r = /7T e (ri = 1/t, t)P(r)dT. By introducing 
the power-law distribution we obtain: 



r + 2r 



t 



7 a 



' T + 2F To 



(3) 



where 7(0:, x) = J* u a ~ 
gamma function. For t 



e u du is the lower incomplete 
00 this expression becomes: 



M«)> T 



r + 2r' 



»(?)' 



(4) 



where T(a) is the gamma function. Therefore, with 
the assumption of a power law distribution for the 
bleaching time r, the fluorescence F(t) emitted by a well- 
defined ensemble of proteins decays asymptotically ac- 
cording to a power law (oc i). As mentioned before, 
this behavior can be understood by considering that at 
time t the term represents the survival probability of 
the system in the fluorescent state. This same argument 
of survival probability allows us to obtain, from the mea- 
sured signal F(t), precious information about the distri- 
bution function of the barriers P{E) without the explicit 
hypothesis of an exponential law. Let us write F(t) in 
terms of an ensemble average of the conditional proba- 
bility P r (r > t\E) to have a lifetime r > t with a barrier 
of height E: 



F(t)=F / P(E)dEP r (r > t\E); 



(5) 



(we note F for F(t — > 0)). The expression J5J can be 
analytically evaluated if we replace the Poisson process 
associated to the Arrhenius law with a deterministic pro- 
cess 0. In this case P r (r > t\E) oc l-®{t-r{E)) with 
noting the Heaviside function, and JSJ becomes 



/>oo />oo 

F(t) =F P(E)dE =F P(E)dE. 

JE(t) Jk B Tln(t/T ) 

This expression can be easily inverted: 

t dF(t) 1 



P(E) 



F dt k B T : 



(6) 



(7) 



where, again, E(t) is such that t = To exp(P/fcsT). We 
can apply such a mathematical transformation to the ex- 
perimental data of Fig. 1 and the result is represented in 
Fig. 3. For energies E > ksT, the distribution deduced 
from the signal (dotted line) is well approximated by an 
exponential law (continuous line). This result strength- 
ens the hypothesis of an exponential distribution of the 
barrier heights which was previously introduced as a pos- 
tulate. For energies E < ksT the retrieved distribution 
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FIG. 3: Barrier height distribution (dotted line) deduced from 
the experimental signal of Fig. 1 by using the expression £J • 
The energy scale on the horizontal axis is deduced in the 
absence of any free parameter. The continuous line represents 
the best exponential fit in the range E > ksT. 



does not follow an exponential law. This behavior is not 
surprising because the Arrhenius law is no more valid 
for t < To (i.e. for E < ksT ). Let us stress that 
the horizontal energy scale in Fig. 3 is obtained in the 
absence of free parameters. Moreover, the order of mag- 
nitude that we obtain for the energy barrier heights is in 
the typical range expected for conformational changes of 
biomolecules. 
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FIG. 4: Parameters a (a) and ta (b) as a function of the 
excitation intensity I. Data are obtained by several experi- 
mental runs on different GFPuv ensembles. While the loga- 
rithmic increase of a is not explained by the theory, the pre- 
factor ta (filled circles in b)) follows an inverse saturation 
law A + B/I as expected from expression 0. The contin- 
uous line represents the best fit to this law (A — 3.7 s and 
B = 0.24 s-MW-cm" 2 ). 



Let us finally analyze the experimental results obtained 
for several excitation intensities I with different ensem- 
bles of proteins. In order to eliminate the possible local 
fluctuations of GFP concentration, we consider the nor- 
malized quantity F n (t) = F(t)/Fo. As previously men- 
tioned, for each I the asymptotic time evolution of F n (t) 
is well described by an empirical law (jA/t) a . The pa- 
rameters ta and a, measured as a function of /, are re- 
ported in Fig. 4. The exponent a increases slightly (with 
a logarithmic law) with I (Fig. 4a). The pre-factor ta 
(filled circles in Fig. 4b) follows an inverse saturation 
law A + B/I (the continuous line in Fig. 4b is the best 
fit to this law). The logarithmic behavior of a is not 



explained by our model in which a only depends on the 
temperature. On the contrary, the expression ac- 
counts for the observed I dependency of the pre-factor 
ta- Indeed, by considering the normalized excited state 
population ;jr7jr^& = r ^ r (•n L e (t)} T , we obtain from ex- 
pression (@J that t2 = a (to Vj ^F ) T{a). By consider- 
ing that r' oc /, we obtain the empirical law that we used 
in Fig. 4b to fit the experimental data. 

In conclusion, we have found experimental evidence of 
the fact that the photobleaching of a well-defined ensem- 
ble of Green Fluorescent Proteins is governed by Levy 
statistics. The model that we have developed is based on 
an irreversible thermally-activated conversion between an 
excited electronic state and a dark state. This model al- 
lows us to explain the observed temporal decay of the flu- 
orescence. Moreover, it gives us quantitative information 
about the probability distribution of the barrier heights 
involved in the aging mechanism. Within this simple 
model the exponent a should only depend on tempera- 
ture. Further experiments, beyond the scope of this pa- 
per, will investigate this T dependence and the observed 
logarithmic / dependence of the a coefficient. 
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